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Abstract 

We study the finite temperature properties of an effective chiral Lagrangian which 
describes nuclear matter. Thermal fluctuations in both the nucleon and the meson 
fields are considered. The logarithmic and square root terms in the effective potential 
are evaluated by expansion and resummation with the result written in terms of the 
exponential integral and the error function, respectively. In the absence of explicit 
chiral symmetry breaking a phase transition restores the symmetry, but when the pion 
has a mass the transition is smooth. The nucleon and meson masses as a functions of 
density and temperature are discussed. 
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1 Introduction 



In previous work |||, |2|, hereinafter referred to as I and II respectively, we have de- 
scribed an effective Lagrangian which incorporates broken scale symmetry in addition 
to spontaneously broken SU{2) chiral symmetry, as suggested by quantum chromody- 
namics (QCD). This Lagrangian contained a potential with logarithmic terms involving 
the glueball field and the chiral a and vr fields. At temperature T = this led to 
a good description of nuclear matter and finite nuclei at the mean field level as well 
as low energy vrA^ scattering data. The extension of this type of model to SU{3) has 
recently been discussed by Papazoglou et al. Q. Here we examine the predictions of 
our SU (2) Lagrangian for the finite temperature, T > 0, properties of hadronic matter 
which are needed in astrophysical applications and in the study of relativistic heavy ion 
collisions. Previous studies |^, ^ of models of this general type at T > have simply 
included temperature effects for the nucleons. Clearly thermal effects for the mesons 
are also significant, particularly those due to the pion which will be dominant at low 
temperatures. 

The analysis at finite temperature is not straightforward due to the logarithmic 
terms in the potential. We have recently suggested in Ref. |^, hereinafter referred to 
as III, a method of treating these terms which involves expansion and resummation 
with the result cast in terms of the exponential integral. This technique was applied 
to the meson sector of the effective Lagrangian. Here we want to carry out a more 
complete calculation by including nucleons as well. As we shall show, this involves a 
square root term which, after analogous treatment, gives results which can be written 
in terms of the error function. The Lagrangian and our thermal analysis are discussed 
in Section 2, with the detailed expressions for the necessary thermal averages being 
relegated to the Appendix. We give our numerical results in Section 3 and Section 4 
contains our conclusions. 
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2 Theory 



2.1 Equations of Motion 

We will simplify the Lagrangian given in II by excluding the isotriplet vector and 
axial vector mesons, the p and ai. They give no mean field contribution to symmetric 
nuclear matter and, since they are relatively heavy, their thermal fiuctuations will not 
play a significant role in chiral symmetry restoration. We can exclude an additional 
term which was introduced to obtain the physical value of the axial coupling constant, 
gA, since it involves the quantity NtN which will not contribute for symmetric nuclear 
matter. We also discard a term e-^NN, which explicitly breaks chiral symmetry, since in 
II tiny values of es in the range to —15 MeV were preferred. An unfavored symmetry 
breaking term labelled e2 is also omitted. Then our effective Lagrangian can be written 
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the explicit symmetry breaking term e'l is present. The field strength tensor for the u) 



field is UJ 



dyUJu- The oj mass is generated by coupling to the glueball field 



as was previously found to be necessary in order to describe nuclei. The mass term 
can be written in terms of the vacuum mass m^^ = G]^^4>q. In Eq. (jl|) we have also 
included a quartic term in the uj field since this generally improved the phenomenology 
in I and II. 

In Eq. (|^) we have replaced the traditional coupling of the pion and the sigma 



meson to the nucleoli, —Ng{a + in ■ t'j^)N, by —Ng^Ja"^ + tv'^N, which is a more 
natural form for including thermal fluctuations. As pointed out by Weinberg such a 
transformation can be achieved by a redefinition of the nucleon field. The form of the 
Lu coupling is invariant, while the extra terms generated by the derivative involve NtN 
and so will not contribute here. (The transformation has no effect if the pion field is 
set to zero as appropriate to T = 0.) 

The quantities B and 6 in Eq. (|l]) are parameters. For the latter, guided by the 
QCD beta function, we take 6 = 4/33 as in previous work. The logarithmic terms 
contribute to the trace anomaly: in addition to the standard contribution from the 
glueball field Q there is a contribution from the a field. There is also a contribution 
from the explicit symmetry breaking, which is related to the pion mass. Specifically 
the trace of the "improved" energy-momentum tensor is = 4evac('?^/</'o)^) where the 
vacuum energy evac = -IB(J)q{1 - 5) - e^. 

We take the vacuum glueball mass to be approximately 1.5 GeV in view of QCD 
sum rule estimates [ p!o| ] of 1.5 GeV and recent lattice estimates |Tl| of 1.7 GeV; small 
shifts in the precise value of the mass are inconsequential here. Since the mass is large 
in comparison to the temperatures of interest, we shall neglect thermal effects for the 
glueball. We define the ratio of the mean field to the vacuum value to be x = 4'/4'o- 
Then the thermal averages of Lagrange's equations for the 0, a and uj fields in uniform 
matter are: 

= 4Box' Inx - Bo6x ^2x^ In (^^) " (^^) ) " ^o<^X^ 
-4X - (^^) ) + e[x' - mlx{u,un , 

= -B,al5x'(^-^^)+ {Bo5x' + e[){a) - e'^^o + (^ J^t-^ NN^, 
= -mlx^{LO^) - AG\ {ojuoj'uj^) +g^{Nj^N), (2) 
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where we have defined Bq = B(j)Q. Here the thermal averages are denoted by angle 
brackets. 

Consider first thermal efi'ects for the a and tt fields. We break a into a mean field 
part a and a fluctuation Aa with mean value (Act) = 0. The mean value of the pion 
field (tt) is, of course, zero. We write 

^ = i^(a-2 + 2a Aa + Aa^ + tt^) 

= 1/2 + 2i^Aiy + , (3) 

where, as in III, u = a/ao, Au = Aa/a^ and ijp = (Ao"^ + 7r2)/crQ. Since expanding 
the fluctuations out to lowest order will not properly treat {a"^ + tt^) when it occurs 
in the denominator or in the logarithm in Eqs. @), we proceed as in III by expanding 
out the cross term 21/ Av of Eq. (Q). The motivation is that at low temperatures the 
thermal average of ■j/'^ is small, while at high temperatures v is small, so the cross term 
is relatively small in both limits. Note that odd powers of Av can be dropped since the 
thermal average gives zero. We also have to contend with products of meson and baryon 
factors. For these we adopt a factorization ansatz, i.e., {f{a, tt)NN) ~ (/(c, 7v)){NN). 
For the omega field we write uj^ = ujq + Auj^ and terms linear in Alj^ will vanish. The 
Gi term, which involves four lo fields, will be treated at mean field level only, since it is 
a small correction and since this avoids the complication of differing longitudinal and 
transverse masses. With these approximations Eqs. (§) become 

+i?ox'(41nx - 5) - e'iX(2z^ - - (^') - x') - rn^xK' + {Au^Au:^')), 
. I 1 2Ai/2 4z/2Ai/2 Sv'^Au'^ \ 
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Z/2 + ^2 (^2+^2)2 (^,2 + ^2)3 (^,2 + ^2)4 / 
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^2 + ^2 (:.2 + ^2^3/2 2(i/2+ ^2^5/2 2(l/2 + ^2^7/2 ^ 



= -mlx^iOQ - AGl^l + 9ujP ■ (4) 

Notice that there is only a solution with exact chiral symmetry, z/ = 0, when there is no 
explicit symmetry breaking, e'l = 0. For most purposes it is sufficient to truncate these 
equations one power lower, however in the absence of explicit symmetry breaking this 
has the disadvantage that there is a small region near the chiral phase transition where 
solutions cannot be obtained. The expansion parameter is of order 4z^^(Ai/^)/(z^^ + 
(i/j^))^ which a posteriori we find to be < 0.05 ~ this is satisfactorily small. 

The evaluation of the thermal average of the logarithm and the integral powers of 
[v"^ + ■0^) was discussed in detail in III. A formal expansion in -0^ was made. The 
thermal average of (i/'^)"' was then written in terms of (■j/'^)" using counting factors 
that assume (V'f ) is independent of i, that is (Ao"^) = (vr^) where tTq is a component of 
the pion field. This is a high temperature approximation which will be accurate when 
chiral symmetry is exactly or approximately restored or when thermal contributions 
are dominant. Nevertheless, our expressions also yield the correct low temperature 
limit. This suggests that our approximation is reasonable, although we do not have 
a quantitative assessment of the errors involved in the intermediate region. The final 
step is to resum the series and the result can be written in terms of the exponential 
integral. A similar procedure is followed for the half-integral powers for which the result 
can be expressed in terms of the error function. This is discussed in Appendix A and 
expressions are listed there for all the thermal averages needed in Eqs. (Q) and the 
equations below. Note that the T = limit of Eqs. @ yields the relativistic mean field 
expressions given in II. 

The evaluation of (i/'^) requires the thermal average of the square of a scalar or 
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pseudoscalar field. The standard results are 



Here P = 1/T is the inverse temperature and iTa represents a component of the pion 



field. The energies, 6,^ = ^/c-^ + m*-^ and Co- = y/k^~+Tn^, depend upon m*^ and m 



*2 



respectively the effective pion and sigma masses. Here the asterisks denote the finite 
temperature and/or density masses which will be calculated in the next subsection. For 
the vector uj the corresponding result is 

oo 


with e^j = \/k'^ + . The nucleon density and the scalar density in Eqs. (^) are 

oo 

p = — J dkk"^ 
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where the effective chemical potential ^* = fj, — QujOJq and the energy E* = VWTm*^ 
with the nucleon effective mass M* defined below. In calculating thermodynamic inte- 
grals, such as these, we find it convenient to make use of the numerical approximation 



scheme of Ref. [12| 



2.2 Masses 

Consistent with our factorization hypothesis, we define the effective nucleon mass to 
be 



M* = g ( Vcr'^ + tt' 



+ ^ -2(.2+ ^2)3/2 -8(,2+^2)V2) ' 
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where we have truncated the expansion at the same order as the equation of motion. 
When r = 0, Eq. (A. 5) indicates that the usual result M* = ga^u is obtained, with 
the vacuum nucleon mass, M = gao, determining g for a given value of ctq. Since 
we consider thermal fluctuations, the nucleon mass will not become zero when chiral 
symmetry is restored and v ^ 0. In fact. 



lim ^ = V „r • (9) 

In the case where e'^ = the sigma and pion arc masslcss at the chiral restoration 
temperature Tc, so (V'^) = T;V(3cr^) and M*/M = 0.543Tc/ao. 

For the meson fields we define the effective mass at finite temperature in terms of the 
thermal average of the second derivative of the Lagrangian. This means that we only 
consider contributions arising from a single interaction vertex. Since the mixing between 
the glueball and the a meson is small, we neglect it here for simplicity. Specifically 



^2/1 
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Notice that at zero temperature a nucleon contribution to the pion mass of goQpsjv is 
automaticahy obtained. In II this required the evaluation of the nucleon loop contribu- 
tion to the pion propagator, whereas here it arises from the form of the coupling used 
for the chiral meson fields and the nucleon. As we have remarked, we do not consider 
thermal fluctuations in the glueball field and so its mass does not enter the equations. 
However it will be useful to display the mass in Sec. 3. The a and tt masses are needed 
in evaluating ("0^) and the mass depends on x- As a result, the equations of motion 
(Q) and the expressions for the masses (10) must be evaluated self-consistently. 



We need to expand the denominators in Eqs. (|10D as discussed in the previous 
subsection. The expressions can be simplified by using the equations of motion (^ for 
the case where 7^ and x 7^ 0; ^ well as the relation 

l/cr\/l \ I \ 

2au{ —-^——) , 11 



do \(cT2+7r2)"/ \(a2+7r2)"/ \ (^2 + 7^2)"+! ^ 

which is valid in our approximation scheme for a an arbitary integer or half-integer. 
We obtain 

2*2 , 4 2/ 1 8Au^ , ^3u^Au'' + 2Au^) \ , e[x' 

\(l/2+'0^)2 (jy2_^^2^3 (j^2_^^2^4 j y 

2 / 1 6Azy2 5(3z^2^^2 ^ 2Av^) \ 



(i,2 + ^2^3/2 (^2 + ^2)5/2 2(l/2 + -02)7/2 



f' v2 
2 *2 ^lA 



4>im%' = 4{Box' + e'li^) - 2{Bod + e[)ii^' + {^p')) + 2^2 (c^2 ^ (Ao^^Au;'^)) . 

(12) 

It is straightforward to verify that in the zero temperature limit the results of II are 
obtained. The explicit symmetry breaking parameter e'l = (cJo^T'Tr)^ is fixed by the 
vacuum pion mass and the chosen value of ctq. For the case when — > the a and vr 
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masses become equal: 



2 *2 2 *2 D e 2 / 1 \ , / Ovr , / 2 no\ 

a^m^ ^a^m^ Bo6x \1 - I + ^^"^•^ U 32(-02^ ' 

When e'l = 0, the expression on the right is vahd for temperatures at which is exactly 
zero. 



2.3 Thermodynamics 

The grand potential per unit volume can easily be written down: 

^ /m 1 *2 2 2 /-(4 4 1 *2/A 2\ 1 *2/ 2\ 

+^ / '^^^'^ [ln(l-e-^^'^)+31n(l-e-^'=")+31n(l-e-^'^^ 
_^ 1 dfc fc2 [in (l + e-^(^*-'^*)) + In (l + e-/3(^*+^*))] . (14) 

The subtraction of the fourth and fifth terms on the right in Eq. ( [I^ ) is necessary to 
avoid double counting |T^ . 

Now if one takes the partial derivative of ^/V with respect to X) or the 
equations of motion (^) are obtained. This is an important and non-trivial consistency 
check. In order to show this one needs the equivalences 

d IJa^ + ^'\\_il 4 \ 



d / V<7^ + 7r" \ a / fJo ^ 

(15) 



\ do / 8 \V(j2 + 7r2, 
We have not succeeded in proving these relations in general, but by expanding out 
2v/S.i' as before and using the explicit equations in the Appendix we have verified Eqs. 
(|l5| ) to orders {v"^ and (i/^-)- ■02^-9/2^ respectively, in our approximation scheme. 
This is all that is needed here. Furthermore Eqs. ([l5| ) are exact in the limits of zero 
and infinite temperature. 
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In practice it is necessary to truncate the expansions. If the equations for the fields, 
the masses and the grand potential are truncated at the same order, then the pion mass 
is exactly zero in the absence of explicit symmetry breaking and there is approximate 
consistency between the equations of motion and the grand potential. Alternatively 
the equations can be truncated at different orders such that there is exact consistency 
between the equations of motion and the grand potential, but then the pion mass is 
non-zero at low temperature when e[ = 0, thus violating Goldstone's theorem. We 
choose the former alternative and note that the inaccuracy is not large, although it can 
result in a small negative value for the pressure at low temperature. 

We therefore truncate the expansion for the thermal average of the potential at 
order {u"^ + i/'^)"^ and write 



We have added a constant term here so that (V) is zero in the vacuum. The pressure 
P is of course —^/V. 

It is straightforward to derive the energy density, which takes the form 



(V) = x'[i?o Inx - \Bo{l + 6) + y,)] + i(5o<5 + e',)x'{i^' + (^')) 




(16) 



/mil *2 2 2 I o/^4 4 1 *2 / A 2\ 
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(17) 
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Table 1: Values of the parameters. 



Quantity 


G4 = 


G^/guj = 0.19 


(MeV) 


236 


228 




10.5 


12.2 


C = (l>o/(^o 


1.28 


1.41 


ao (MeV) 


110 


102 


e;^/^ (MeV) 


123 


119 



3 Results 

The parameters used to obtain the numerical results are listed in Table 1. They were 
determined in II by fitting to the properties of nuclear matter and finite nuclei. Quite 
a reasonable phenomenology was obtained there, with the non-zero value of G4 being 
slightly favored, so it is sensible to explore the behavior at T > 0. We shall also consider 
the case where the vacTium pion mass vanishes (e'^ = 0) for which the parameters differ 
only slightly from those specified in Table 1. 

If e'l = 0, so that the pion is massless in the vacuum, chiral symmetry can be ex- 
actly restored at sufficiently high temperature. We show in Fig. 1 the chiral restoration 
temperature Tc versus the ratio of baryon density p to the density at saturation pQ. The 
primary distinction between the zero and non-zero G4 cases arises from the different 
values of ctq in Table 1, which fix the overall temperature scale. Thus the same qual- 
itative behavior is seen in the two cases with the corresponding temperatures slightly 
larger when G4 = 0; at zero density this scaling factor is almost exactly 110/102 ~ 1.1, 
but it increases to about 1.7 at very high density. Therefore it is sufficient in the 
following to focus on the G4 ^ results. 
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In III, baryons and uj mesons were absent and Tc was 187 MeV. Since at restoration 
m* = m* = 0, Eq. ([l3|) shows that Tc is modified by the scalar density which is non- 
zero even at zero baryon density. We obtain a value of 162 MeV, which is a reasonable 
order of magnitude in view of the many estimates in the literature. We mention that 
the standard Gell-Mann-Levy model |l|] yields Tc = V^fn = 132 MeV. Also Gerber 



and Leutwyler |15| estimate Tc = 190 MeV in two- flavor QCD using an effective chiral 
Lagrangian to three loop order, while a recent lattice calculation |l6[ gave 140 — 150 
MeV. The restoration temperature decreases with increasing density so that at po it is 
125 MeV and asymptotically it is 45 MeV. A finite Tc is required in the high density 
limit since at T = 0, as shown in Fig. 2, z/ becomes ~ 0.3, which is small but non- 
zero. Thus chiral symmetry is not completely restored at zero temperature in our 
model; this finds some support from lattice calculations 0]. Papazoglou et al. |^ 
have found that there are regions in parameter space for models of this type where 
chiral symmetry is restored at high density, but with these potentials the compression 
modulus is unreasonably high. 

The remainder of Fig. 2 shows that u becomes zero at progressively lower densities 
as the temperature is raised. The dotted lines here indicate metastable regions and, 
while the transition is clearly first order, we do not believe that our approximation 
scheme is sufficiently accurate to reliably determine the order. For instance, if our 
expansions are truncated at a lower order we are unable to find solutions for a small 
region in the vicinity of the transition. Suppressing results in the metastable regions, 
we plot the sigma mass in Fig. 3. For T = the mass increases with density due to 



the scalar density contribution in Eq. (12). With increasing temperature the familiar 
behavior is seen, namely, the mass drops to zero to become degenerate with the pion 
mass at successively lower densities. Beyond this point the two masses increase together 
as either the temperature or density is raised. 
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We now turn to the more realistic case where ('^ ^ and the pion has its physical 
mass in the vacuum. Figure 4 plots v as a, function of density. At low temperatures 
V increases slightly for finite densities, before decreasing as the temperature is further 
raised and becoming rather independent of density. Substantial reduction in v requires 
higher temperatures than needed for Fig. 2, which seems reasonable in view of the 
finite pion mass here. With baryons present v smoothly tends towards zero as T — > oo, 
whereas in their absence it stabilized at ~0.15, as demonstrated in III. In either case 
there is no phase transition when e'l ^ 0. In this connection it is worth mentioning 



that the lattice QCD calculations of Brown et al. [18| show that a second order phase 
transition for massless quarks is washed out when they are given a finite mass. 

The pion and sigma masses are shown in Figs. 5 and 6. As the temperature is 
raised the pion mass increases, while the sigma mass decreases so that they approach 
each other. They are degenerate when a temperature of ~ 250 MeV is reached and 
depend very little on density. Thereafter they continue to rise, passing the 1 GeV mark 
at T ~ 300 MeV (not illustrated). The sigma mass remains above 400 MeV for all 
densities and temperatures. This is qualitatively different from the high temperature 
behavior for the e'^ = case shown in Fig. 3. Furthermore, comparison with III 
indicates that the presence of baryons pushes up the mass noticeably. 

The nucleon effective mass, M*, in shown in Fig. 7. At low temperatures its be- 
havior is dominated by v and the mass increases slightly with temperature for non-zero 
densities. However, the thermal fluctuations, (i/'^), begin to dominate for temperatures 
above 150 MeV, leading to a drop in M* at low densities. This thermal domination 
implies that the behavior of the mean field v is of minor consequence, which is why a 
plot of M* for the e'l = case looks quite similar to Fig. 7. By T = 250 MeV, M* 
is rather independent of density. The mass continues to fall at higher temperatures 
and remains roughly independent of density. At such high T this is not caused by the 
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very small but rather by the diminished value of {ijj^) which follows from the sharp 
increase in the meson masses seen in Figs. 5 and 6. Qualitatively the overall behavior 
of M* is strikingly similar to that reported by Furnstahl and Serot for the Walecka 
model. For a given T, their mass drops to a lower level than ours, probably because 
they did not consider meson fluctuations. On the other hand at least one model - the 
quark meson coupling model - shows a different behavior |]20t| in that the nucleon mass 
simply increases monotonically with temperature. 

The ratio of the glueball mean field to the vacuum value, x = 4'/4'Q: is shown in Fig. 
8. This ratio differs from unity by more than 10% only for the highest temperature, 
250 MeV, and low to moderate densities. To complete the discussion, we must examine 
the glueball effective mass of Eq. (^), m^, which is displayed in Fig. 9. Again we 
remark that this mass was not used in the self-consistent solution of the equations 
since glueball fluctuations were excluded on the grounds that the mass is large. This 
is proven justified for low temperatures, but it is becoming questionable at T = 250 
MeV for low density. Also the difference between and m* starts to become small 
and it may not be adequate to ignore the mixing induced by the off-diagonal terms 
in the mass matrix described in II (the vacuum mass in Fig. 9 will be pushed up a 
little by this mixing). Furthermore, beyond this temperature the fluctuations of the uj 
meson field, (Auj^iAuj^), start to become significant. This reduces the glueball effective 
mass further and continuing to solve for m* and x self-consistently quickly sends 
to zero at T ~ 300 MeV. This effect was not observed in III since the uj meson was 
not included. To obtain reliable solutions in this region one would first of all need 
to consider fluctuations in the glueball field and compute the mixing with the a field. 
Secondly, since the mass of the p meson is similar to that of the uj, the fluctuations 
in this fleld will also play a role even though its mean fleld is zero. In addition one 
should include the chiral partner of the p, namely the ai, since with the approximate 
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restoration of chiral symmetry the masses will be similar. If the effect of fluctuations in 
the p and ai fields is estimated by appropriately increasing the degeneracy factor for the 
cj fluctuations, it is found that m*^ is driven to zero at a temperature of ~200 MeV (at 
zero density; the effect is softened with increasing density, but still present). Although 
we expect that glueball fluctuations will partially counter this effect, it is intriguing 
since it could be an intimation of deconfinement. Lattice calculations suggest that this 
occurs at a similar temperature to chiral restoration [^]. However a consistent thermal 
analysis for the (j)^ p and ai mesons would be quite involved and is beyond the scope of 
the present work. 

Turning to the thermodynamics of our model, we first remark that in common 
with other models we observe a liquid-gas phase transition; the critical temperature is 
found to be 16 MeV. The pressure and energy density are plotted in Figs. 10 and 11, 
respectively, as functions of density. The calculations here are for e'^ ^ 0, but taking 
e'l = yields similar results, except that all numbers are reduced by a few percent in 
the absence of a finite pion mass. At the highest temperatures and densities shown, 
it can be noted that E ~ 3Py, which is the massless gas result. However this does 
not hold for lower values of either T or /), for which interactions dominate the thermal 
fluctuations. For orientation we have compared the pressure to that obtained with a 
crude lowest order treatment of gluons and massless quarks, taking the bag constant to 
be |evac|- We find that at no point, either in temperature or density, is the quark-gluon 
plasma the preferred phase prior to chiral restoration (for e'l =0). At zero density, 
for example, chiral restoration is achieved at Tc = 162 MeV while the quark-gluon 
plasma has a larger pressure for T ~ 170 MeV. At equilibrium density, where we found 
Tc = 125 MeV, the quark-gluon phase is preferred at T = 164 MeV. 
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4 Conclusions 



We have discussed the finite temperature behavior of an effective Lagrangian with 
whicli we have successfully described nuclear matter and finite nuclei at T = 0. With 
nucleons present, we have thermally averaged a square root term involving the meson 
fields in addition to a logarithmic term. The latter was handled in III by expansion and 
resummation of an infinite series with the final result cast in terms of the exponential 
integral. For the former the result of a similar approach was written in terms of the 
error function. 

Our results showed that at sufficiently high temperature the mean value of the 
a field became small, signalling chiral restoration. In the absence of explicit chiral 
symmetry breaking, a phase transition restored the symmetry at temperatures ranging 
from 162 MeV at zero density to 45 MeV at very large density. It was estimated that 
that the quark-gluon phase would be preferred at somewhat higher temperatures than 
these. In the physical case, where explicit chiral symmetry breaking was present and 
the pion had a vacuum mass, a smooth restoration of chiral symmetry was found. The 
masses of the pion and a meson were virtually degenerate and the order parameter i/ 
became small for temperatures of ~ 250 MeV. At such a temperature the nucleon mass 
was also reduced, but the effect of meson fluctuations yielded some stabilization so that 
M* /M ~ 0.43. At a similar temperature the glueball mass started to be significantly 
reduced due to the coupling to lo vector meson fluctuations. This could be a hint of 
interesting physics, but we were unable to track it further with the present formalism 
since glueball fluctuations and the p and ai mesons are expected to be of significance. 

We thank S. Rudaz for useful discussions. We acknowledge partial support from 
the Department of Energy under grant No. DE-FG02-87ER40328. G.W.C. thanks the 
University of Minnesota for a Doctoral Dissertation Fellowship. A grant for computing 
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Appendix A. Thermal Averages 



A.l. Half-Integral Powers 

We first write the expansion 



1^2 + ^2 



1^(1 + 



2u^ 



E(-i 



2n ' 



n=2 



,n (2n-3)!! ^ 
2"n! i/2n 



(A.l) 



For the purposes of evaluating the counting, we assume that the thermal average (V'f ) 
is independent of the label i, which amounts to assuming the masses of the particles 
involved are the same. Then the result of taking the thermal average of each possible 
pair of fields at a general vertex can be written ((i/?^)") = c„(?/'^)" and, as shown in 
III, c„ = (n + l)!/2". Using this result and defining 
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we have 
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This is a divergent series which we regard as a formal expansion and it must be re- 
summed before it can be evaluated. Matching this expression to the asymptotic expan- 
sions of i"erfc(2), the repeated integrals of the complementary error function [^], we 
have 



Zy2 + -02 



42 



3iv 
4^ 



2zi^erfc(z) + 3i°erfc(z) 



2z + {l 



2 „2 



I— 2 

)\/7re^ erfc(z) 



Here erfc(z) = 1 — 27r 2 e * dt. In the limit of low temperature {z 00) 



V2 + 0^ 



2u 



(A.4) 



(A.5) 
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as follows from expanding out the square root directly. In the high temperature limit 



3 /F 
4V2 



(A.6) 



and the numerical factor is 0.9400 which is close to unity, as one might expect. 

The other expressions that are needed can be obtained by differentiating with re- 
spect to v"^ and by using the relations 



(A.7) 

which are easily shown for an arbitrary function /, when is independent of i. One 
obtains the following: 



+ -02)-^ \ = |_ J2z + (1 - 2z2)^e^'erfc(^) 



(i.2 + ^2^-i\ = _f_ [2z - (1 + 2z2)0Fe^'erfc(z) 



Az/2(i/2 + ^2)-f \±l2z + Az^ + (1 - - 4/) V^e^'erfc(z) 
/ Si' L 



Az.2(^2^^2^~| 



Z 



4 



lOz + 4z^ - (3 + 12^2 + 4z'')V7re'= erfc(z) 



30z/5 



+ ISz^ + 4/ 



1201^'^ 



-(15z + 20z^ + 4/)V7re'='erfc(z) 
66z + 56z^ + 82;^ 

-(15 + 90^2 + eOz"^ + 8z^) V7re'='erfc(z)l . (A.^ 



A. 2. Integral Powers 

The thermal averages needed can be derived in the manner indicated above using 
the thermal average of the logarithm. This is evaluated in analogous fashion to the 
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square root and the details were discussed in III. We simply list the results here, in 
terms of the variable y = z^: 

+ ^2)^ = + (1 _ y)eyEi{y) + 1 , 

y 



r 

,2 



[l-(l + y)e^i?i(y)] 



(^2 ^ ^2)-3^ = A_ ^1 + y _ (2y + y^)eyE^{y) 



Az.2(^2^^2)-4 



4zy2 



24i/6 



l + y-(2y + y2)eS/^,(y) 
3 + y-(2 + 4y + y2)e2'^i(y) 
2 + 5y + - (6y + 6y^ + /)e^^i(y) 



11 + 8y + y^ 



-(6 + 18y + V + /)e^^i(y) . 



The exponential integral is defined |2^ by Ei{y) = dtt 



(A.9) 
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175 




Figure 1. The chiral restoration temperature Tc as a function of the ratio of the 
density to equilibrium nuclear matter density, p/po, for two values of G4 with e[ = 0. 
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Figure 2. The mean sigma field, v = a/ao, as a function of density for various 
temperatures (in MeV) with e'l = 0. The dotted curves indicate thermodynamically 
unstable regions. 
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Figure 3. The sigma effective mass, m* , as a function of density for various tempera- 
tures (in MeV) with e'l = 0. 
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Figure 5. The pion effective mass, m*, as a function of density for various tempera- 
tures (in MeV) with e[ > 0. 
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Figure 7. The nucleon effective mass, M*, as a function of density for various tem- 
peratures (in MeV) with e[ > 0. 
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Figure 8. The glueball mean field, % 
temperatures (in MeV) with > 0. 
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Figure 9. The glueball effective mass, m^, as a function of density for various tem- 
peratures (in MeV) with e'l > 0. 
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Figure 10. The pressure as a function of density for various temperatures (in MeV) 
with e'l > 0. 
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Figure 1 1 . The energy density as a function of baryon density for various temperatures 
(in MeV) witli e\ > 0. 
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